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\Q | Abstract. The if— matrix, also known as the "Wigner reaction matrix" in nuclear 

scattering or "impedance matrix" in the electromagnetic wave scattering, is given 
essentially by an M x M diagonal block of the resolvent (E — iJ) _1 of a Hamiltonian 
i-C | H. For chaotic quantum systems the Hamiltonian H can be modelled by random 

Hcrmitian N x N matrices taken from invariant ensembles with the Dyson symmetry 
i-G , index j3 = 1, 2 or 4. For (3 — 2 we prove by explicit calculation a universality conjecture 

& ' by P. Brouwer [13] which is equivalent to the claim that the probability distribution of 



K, for a broad class of invariant ensembles of random Hermitian matrices H , converges 
to a matrix Cauchy distribution with density ViK ) ex [det (A 2 + (K — e) 2 )] in the 
limit N — > oo, provided the parameter M is fixed and the spectral parameter E is 
^- . taken within the support of the eigenvalue distribution of H. In particular, we show 

that for a broad class of unitary invariant ensembles of random matrices finite diagonal 
blocks of the resolvent are Cauchy distributed. The cases f3 = 1 and {3 — 4 remain 



t^J- , outstanding. 
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The phenomenon of chaotic resonance scattering of quantum waves (or their 
classical analogues) has attracted considerable theoretical and experimental interest for 
more than two decades, see e.g. articles in pQ. The resonances manifest themselves via 
fluctuating structures in scattering or transport observables, and statistical properties 
of such objects can be successfully described within the framework of the Random 
Matrix Theory [2], [3]. The most important object in such an approach is the energy- 
dependent M x M random unitary scattering matrix S(E), S Jf (E)S(E) = ljv/ which 
relates amplitudes of incoming and outgoing waves. Here the integer M stands for the 
number of open channels at given energy, the dagger denotes the Hermitian conjugation 
and 1m is the M x M identity matrix. Statistical properties of scattering observables 
considered at a fixed energy E of incoming waves can be inferred from the corresponding 
probability density of S = S(E) derived starting from rather general physical principles. 
Those include unitarity, causality and (if relevant) the time-reversal invariance imposed 
on S combined with the assumption of maximal entropy (minimum information). The 
procedure yields the so-called Poisson's kernel distribution wih density pi]: 



G /3 






det[l M -S'S] 



det[l M - S ] S] 



-in, 2 



Q3M+2-/3)/2 



(1) 



where S stands for the mean of the scattering matrix, /3 = 1,2,4 is the parameter 
related to underlying symmetries with respect to time reversal and Cp is a normalization 
constant. The mean S is determined by the details of coupling of the systems to 
continuum and thus contains all information which should be specified for a given 
scattering system. In particular, for the simplest, yet most fundamentally important 
"perfect coupling" case S = 0, and the density in (CQ) is constant, implying that the 
S— matrix is uniformly distributed over the unitary matrices of given symmetry. 

Although the above method has proved to be very successful in the statistical 
description of scattering characteristics at fixed energy [2] , it can not be used to study 
statistics of fluctuations of the scattering observables over an energy interval comparable 
with a typical separation between resonances. The latter task can be most successfully 
achieved in an alternative powerful approach going back to the pioneering work [5] which 
is based on the paradigm of random matrix properties of the underlying Hamiltonian 
H describing quantum chaotic behaviour of the closed counterpart of the scattering 
system. In such an approach the resonance part of the S-matrix is expressed in terms 
of the resolvent of such a Hamiltonian as 

S(E) = \ M ~ iff , K{E) = W\E-H)^W, (2) 

where W is an NxM matrix of energy-independent coupling amplitudes between N 
energy levels of the closed system and M open scattering channels. To study quantum 
chaos-induced fluctuations of S one then replaces the Hamiltonian H with a random 
matrix taken from one of the standard random matrix ensembles, usually Gaussian 
Unitary (GUE, — 2) if one is interested in the systems with broken time reversal 
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invariance or Gaussian orthogonal (GOE, (3 = 1) if such invariance is preserved, the 
case (3 = 4 being relevant for systems with spin-orbit scattering. The approach proved 
to be extremely successful, and quite a few scattering characteristics were thoroughly 
investigated in that framework in the last two decades, mainly by the supersymmetry 
method [5j- [Tj . The results of such calculations are found in general to be in good 
agreement with available experiments in chaotic electromagnetic resonators ( "microwave 
billiards") or acoustic reverberation cameras, see e.g. [SJ [9], [10] and most recently in 
[TT] , as well as with numerical simulations of scattering in such paradigmatic model as 
quantum chaotic graphs [T2] . 

The two random matrix approaches described above look very different in their 
formulation, yet they are meant to describe precisely the same object, the S— matrix for 
a chaotic system. The consistency therefore requires that the Poisson kernel distribution 
(PQ) for S must follow from the law of distribution of H entering the relation (T5]). 
Surprisingly, a direct verification of such a correspondence turns out to be a rather 
challenging task. The challenge here is that the two objects are related via the resolvent- 
like K— matrix, and to convert the law of distribution of H into that of the resolvent 
is not at all trivial. A very elegant indirect way round this problem was discovered 
by P. Brouwer [13] who proposed to choose H from the Cauchy ensemble of random 
matrices with density P(H) oc det [A 2 + (H — e) 2 ] , where A, e are two real 

parameters. The main advantage of such a choice is that the resolvent (E — H)~ l is 
Cauchy-distributed as well albeit with modified parameters and, moreover, diagonal 
blocks of Cauchy matrices have again closely related distributions. Using these facts 
Brouwer indeed was able to demonstrate the validity of the Poisson kernel for such a 
choice of H for all values of f3 = 1, 2, 4. He then showed that in the large- N limit the 
eigenvalue correlation functions in the Cauchy ensemble (called "Lorentzian" ensemble 
by Brouwer) have the standard Dyson form, and conjectured that such equivalence 
of eigenvalue correlation functions should be enough to ensure the same S— matrix 
distribution is to be shared by all representatives of the corresponding universality class. 
Although such conjecture sounds very natural, the particular mechanism by which the 
generic spectral properties of H are translated into universality of the probability density 
of the K— matrix and then P(S) remained unclear. To the best of our knowledge no 
further attempts to verify universality of the S— matrix distribution were undertaken 
in the literature apart from (i) the simplest case M = 1 and H G GUE considered 
in [7] and (ii) the recent work [H] which however concentrated on the universality of 
two-point spectral correlations of the individual S— matrix entries rather than on the 
one-point matrix distribution. 

In the present work we establish the universality of the Poisson distribution for 
the S'-matrix under the condition of equivalent coupling to continuum in all scattering 
channels. Since the K- and .S-matrices are related via the Cayley transformation fl2]), 
the claim of the Poisson kernel distribution for the S— matrix is equivalent to claiming 
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that the K-matrix is Cauchy with density 

Pp(K) oc [det (A 2 + (K- e) 2 )] ~^ M+2 -V f (3) 

where the width A and the mean e of the Cauchy distribution are determined by the 
mean S of the S'-matrix and vice versa. Such equivalence can be verified by transforming 
to the eigenvalues of the S'-and f^-matrices. 

In the present work, under fairly generic assumptions on H, we verify that the law 
of distribution of the i^-matrix is indeed Cauchy and relate its parameters A and e to 
the strength of the coupling amplitudes W and the density of states of the underlying 
matrix H as well as details of its (invariant) distribution. We would like to note 
that the statistical characteristics of the i^-matrix are directly accessible in microwave 
experiments [9J where it is related to the real part of the impedance in the regime 
of small losses. Our paper stems from attempts to understand better the mechanism 
behind universality of the probability distribution of finite blocks of random matrix 
resolvents and to provide an ab initio explicit derivation of this distribution for generic 
invariant ensembles of random matrices. 

We will start with showing that the universality of the fC-matrix in question follows 
from the universal limit of a very general spectral object - the product of the ratios of 
powers of characteristic polynomials det(E — H) of random matrices H. 

To that end it is necessary to mention that in the literature there exist two 
alternative choices of the matrix W of coupling amplitudes in (J2J). The standard choice 
is to follow the original paper [5] and to consider the columns Wi,...Wm of W as 
mutually orthogonal iV— component vectors, real for /3 = 1 and complex for (3 = 2. The 
case of equivalent channels then corresponds to (w a , w fe ) = jS a b for all a, b = 1, . . . M, 
or, equivalent ly, 

W^W = 7 1m, (4) 

where 7 > is a coupling constant. An alternative choice which was suggested originally 
in [15] is to consider the columns of W to be independent Gaussian vectors with joint 
probability density function 

P(iy)«e-f wtlf , {W*W) = 7 1 M . (5) 

Both choices are expected to lead to the same results in the limit N — > 00 as long as 
the number of channels M remains fixed. Such an equivalence was explicitly verified 
in [16] for particular scattering characteristics (Wigner delay times), but is expected to 
hold generally. 

We shall first consider the random amplitude case (jSJ) and show the equivalence 
to the fixed amplitude case (HI) for (3 = 2 at the end of the paper. For notational 
convenience, it is more convenient to work with the rescaled ii'-matrix 

K=-^— K = W\E-H)- 1 W, 
7rp{E) 



Universal K— matrix distribution in (3 = 2 Ensembles of Random Matrices 5 

rather than the i^-matrix itself, with p(E) being the large- iV limit of the mean eigenvalue 
density of H at point E inside the support of p(E) (so that p(E) > 0). 

We start with the representation P(K) = f J-p^{X) exp(i|| TtXK) dX of the 
probability density function of K in terms of the characteristic function 



jFj3,n{X) = ( exp 



2 J 



(6) 



where the matrix X has the same dimensions and symmetries as the matrix K and 
the angular brackets stand for the averaging over all random variables the matrix K 
depends on. Writing K = - h^ W^RW, R = (E — H)' 1 , we first perform the averaging 
over the coupling matrix W in <^j which amounts to performing a Gaussian integral 
over the vectors w c : 



-4(-TrWtW+ 



Ttp(E) 



TrXW^RW) 



M 



dW oc W 



c=l 



det 1 N + 



l^Xc 



nNp(E) 



R 



(7) 



Here x\,...,xm are the eigenvalues of the Hermitian M x M matrix X and dW 
stands for the appropriately normalized Lebesgue measure on the space of complex 
or real N x M matrices W. The easiest way to verify ([Z|) is by diagonalizing 
X = T diag(xi, . . . ,Xm) T^ 1 where T is orthogonal for j5 = 1 and unitary for (3 = 2. 
Then one changes (TW) — > W and exploits the invariance of W^W and the measure 
dW with respect of such a transformation. At the next step we bring the characteristic 
function Fp tN (X) to the following form: 



M 



^, N (x) = ( n 



[det (E - H) 



i [det (E + - 



i-,.i r 



Np(E) 



-H 




(8) 



where the angular brackets now stand for the averaging over the N x N matrices H. 
The above relation is exact in the random amplitude model (jSJ) for any choice of iV 
and M. We will show at the end of the paper that for (3 = 2 the same equation (|Sj) is 
valid asymptotically in the fixed amplitude model (Hj) in the limit N ^> M provided the 
probability density of H is rotationally invariant. 

For (3 = 2 and (3 = 4 the object in the right-hand side of (jHJ) is well-studied in the 
Random Matrix Theory [TT1 [TBI fT9] . In particular, in the simplest case (3 = 2 the formula 
(2.14) from [18j appears to be most useful for our goals. Namely, for N x N matrices 
H distributed according to an invariant ensemble density with polynomial potential V, 



P(H) oc exp [-NTrV(H)], V(H) = ^Q# 2 ', 



c p > 0, 



(9) 



1=0 



(10) 



" 1)/2 ex P ( ~ na E $^(Cc - Vc) j A2 | C | A2 ^| det(S(Q - Vj)), 
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the following universal relation holds asymptoticalljl||: 

/^ det(E + r]c /(Np(E))-H) \ 
n^oo \l\ det(E + ( c /(Np(E)) - H) ] h 

^M(M-l)/2 exp j _ 7ra£ ^ _ ,\ A{C, V} 

where A{rj} = Ili<i(^i — Vj) * s the Vandermonde determinant, and 

ctf „\ _ exp(i7rsgn(ImC)(C-^)) n _ V\E) 
C - V 2Tvp{E) 

An analogous result for averaged products of ratios of characteristic polynomials with 
(3 = 4 is also known [19], but has a more complex structure, with Pfafnans replacing 
determinants. Unfortunately, for /3 = 1 no result of comparable generality seems to be 
known for the products of square roots of the characteristic polynomials , though for 
M = 1, 2 it can be in fact evaluated in closed form, see e.g. [21] and references therein. 
Below we consider in full generality only the case of Hermitian ensembles with (3 = 2, 
whereas the cases (3 = 4 and especially (3 = 1 remain a challenge to us and are currently 
under investigation. 

With the asymptotic relation (fTUj) in hand, one can evaluate the characteristic 
function ([8]) of the rescaled A-matrix in the limit iV — > oo and M fixed. 

Proposition 1 Assume that the N x N matrix H has invariant distribution (TJ|). Then 
in the random amplitude model (TJJ) we have linijv-s.oo*7 : /3=2,7v(A) = J r g =2 (X) ; where 

(gM-i{xi) ... 9m-i(xm)\ 



T^ 2 (X) = H M(M-m eM-p^TrX] ^ 

A{X\ 



(12) 



\ 9o{xi) ... g (x M ) J 



with 

M-n 



g M -n(x) =exp(-7|x|)a; n x ^ f[W- ( 13 ) 

1=0 

Proof. To adjust equation (fTUj) to our goals we first set r/i = r] 2 = . . . = r\ M = there. 
In this limit A{(,r]}/(A{(}A{r]}) — > (d x . . . x Cm) M and equation ([TO]) becomes 

tL (n det( d ;f^' h)) H = exp {t ^ i:sgn(lm « - -«) 

(r r , M 1 [g(d,vi) ■■■ 9(Ci,vm)\ 

x ( \M(M-l)/2 Isi x • • • x Cm) .. I 



, r , l im a r i det 

A{C} »?i...i7M->o A{r/} 



VKCm,*7i) ••• 9{Cm,Vm)J 



(14) 



| We restrict ourselves to the polynomial potentials in (j9|) for the notational convenience. The 
asymptotic relation (|10|) . and as a consequence our Proposition[T]hold for invariant ensembles of random 
matrices under fairly general conditions on the matrix measure, see the recent paper [20 . 
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where 



g{(,v) 



exp[— m sgn(Im Qrj\ 



(15) 



The limits are now performed successively applying the L'Hospital's rule, the final result 
for the second line in equation ( Tbf|) being 



lim - 



det[sf(Ci,%)] 



l<i,j<M 



,n=l 



fifo(Ci) \ 



(16) 



\<7m-i(Cm) ••• 9o((m)J 



where we have defined 



~9n(0 



d" 



drj'' 



-g((,v) 



ni 



^ S(-"0! 



[— z7rsgn(ImC)] n C 



n-l/—l-\ 



;i7) 



Finally, by redefining g n {C) = e ln ^ ssn( - lm -^( M g n (()/n\, several factors in front of the 
determinant can be absorbed into the determinant. After identifying ( c — > i^x c /ti we 
arrive at (fT2"|) . This completes our proof of Proposition [TJ 

At the next step we observe that achieving our main goal is equivalent to verifying 
that J r /3 = 2(X) is the characteristic function of a matrix Cauchy distribution. 



Proposition 2 



-i Tr KX 



dK 



oc ^ =2 (A) 



det [ 7 2 + (K- 7 a s ) 2 ] M 
where the integral is over the set of all Hermitian M x M matrices K . 



(18) 



Proof. A standard random matrix calculation which involves changing the variables of 
integration in f )18p to the eigenvalues and eigenvectors of K and then applying the 
Itzykson-Zuber-Harish-Chandra (IZHC) formula, see e.g. [22J, and the Andreief-de 
Bruijn integration formula yields 



/ 



-iTrKX 



dK 



-i^/ciE TrX 



det[ 7 2 + (K - 7 a B ) 



21 M 



OC 



A{X} 



det 



f'(x 1 ) 
Kf^-'Kxi) 



f(x M ) \ 

f (M - l Kx M )) 



(19) 



Here 



/(*) 



-ikx 



dk 



( 7 2 + k 2 ) M 



ex \>yz\ M -*K M _i{\>yz\), f (M \x) 



Mi 



iM 



dx M 



M 



(20) 



where K v {x) is the modified Bessel (Macdonald) function. In particular, for M = 1 we 
have f(x) oc e^ 7 ^', for higher M we have 



frM^g iijg^fe M 



M-i-; 



(21) 



Universal K— matrix distribution in j3 = 2 Ensembles of Random Matrices 



where we added a subscript to indicate the M- dependence. Using a recursive relation 



for the derivative of the Macdonald function we can show that f' A 
Inductively one gets for higher derivatives 



M 



X 



I 



M ' 



x) 



\m/2\ 

£ 

1=0 



m!(-l) 



m—l 



l\(m-2l)\2 l 



(lx) m - 2l f M - m+ i(x) 



-lxf M ~i{x). 



(22) 



where ['J denotes the floor-function. This enables us to simplify the determinantal 
structure of ( TT9l) . By successively adding to the n-th row appropriate linear 
combinations of all preceding rows 1, 2, . . . , n — 1, and exploiting yet another recursive 
relation for the Macdonald function one can remove all terms in the equation (J22l) but 
the one for / = 0, leading to 



det 



/ /w(a?i) 
/mOi) 



\f?-%£ 



Im(xm) 



\ 



f^MJ 



oc det 



( ImM 

xif M -i{xi) 



!m{xm) \ 

XmIm-I^Xm) 



(23) 



Vxf-Vi(xi) ... xZ^hM) 



where the combination involved in the n-th row in the right-hand side is given explicitly 
by 

{2M-2n-l)\ 



„n-l 



./-' f M -n+i(x) oc e ' -'./••■ 



-7|x| 



M-n 
1=0 



U{M-n-l)\2 M - n - 1 



\ix\ 



(24) 



The equations (TTBl and (1241 have a very similar structure, though the coefficients of 
the terms in the sum are still different. In fact this similarity can be further exploited 
to show that the determinants in equations ffT2l) and (fT9j) (or equivalently (T23]) ) are 
proportional to each other, thus verifying the equation ( TT8|) . 

We start our demonstration of this fact with bringing the first row of the 
determinant in equation ( TT2l) to the form coinciding with the first row of the determinant 
in ( 1231) . Since the zeroth and the first order coefficients of gM-i(x) are both equal to 
unity, and the two corresponding coefficients are also equal in the expression for /m(x) 
(but are different from unity) we can safely change those coefficients in gM~i{x) to the 
coefficients in Jm{x) as such a change gives rise to a constant proportionality factor for 
the determinant. 

The main observation is that the adjustment of both the coefficients a n and a n +i, 
given that all previous coefficients are already adjusted, can be done simultaneously by 
adding the (2n + 1)— th row multiplied with the factor 



-1^ 



(2M -2n-2) 



(25) 



n\(M -n-l)\2 M - 1 ' 
For this procedure to work we need to verify, for any integer n, the following identity 



1=0 



(2M-2/-2)! 



(2M-2n- 6-2)1 



l\(M -I- 1)!2 M ~ 1 (2ra + 5 - 2l)\ (2ra + 8)\{M - 2n - \)\2M-^-s-i 



(26) 
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with 5 = or 5 = 1. The left-hand side of (1261) is what becomes of the (n + 5)-th order 
coefficient of gM-ife) after adding multiples of all odd rows up to 2n + 1, choosing the 
multiplication factors according to ( 1231) . The right-hand side equals to the corresponding 
(n + 5)-th order coefficient of Jm{%)- Both equations can be conveniently combined into 
a single relation: 

where M. = M — 1 and m = 2n or m = 2n + 1. To verify (1271) we first express the first 
binomial on the left-hand side by a contour integral using its generating function and 
the Cauchy's residue theorem. The summation over I is then performed in the integrand 
using the binomial theorem, and the resulting contour integral can be again evaluated 
by the residues, yielding precisely the right-hand side of the relation (l2"Tj) . 

We conclude that it is indeed possible to transform gM-i{x) into /m(^) by adding 
multiples of all odd rows to the first row. Note that in each step two coefficients 
get adjusted simultaneously, and this is precisely the mechanism ensuring the whole 
procedure being functional. Had it not been for that property, we would be only able to 
change half of the coefficients to the required form, since adding even rows to odd rows 
or vice versa is meaningless due to their rather different structure. All remaining odd 
rows as well as all even rows can be treated by exactly the same procedure, since the 
coefficients involved are essentially the same as before. Note also that as the very last row 
contains on both sides the function e~ r/ ^x n ~ 1 the coincidence is ensured automatically. 
This completes our proof of Proposition [2l 

Since the characteristic function uniquely determines the law of distribution, one 
concludes from Propositions [Hand [2] that the distribution of the i^-matrix (J2J) converges 
in the limit N — > oo to the matrix Cauchy distribution with density Pp = 2{K) (jHJ) having 
mean e = r yV'(E)/2 and width A = n r yp(E). This corresponds to the Poisson kernel 
distribution (TJQ) for the S-matrix with mean Sij = 1 7^ 7P L][ 1 l^ g ] ^j- The case of perfect 
coupling is then obtained for a(E max ) = and n^p(E max ) = 1, where E max denotes the 
point where p(E) has its maximum. Thus indeed, the Poisson kernel distribution for 
the S'-matrix is universal in the random amplitude model (JHJ) in that it does not depend 
on the choice of the random matrix ensemble for the underlying matrix H. 

Finally, we would like to demonstrate that the fixed amplitude model fllj) yields 
the same universal behaviour of the K— matrix in the limit N — > oo. Let us again 



consider the characteristic function J z j3=2,n{X) = ^exp[— iTr( ^l^ XW^RW)] 

(exp(-iTiT x UR A W)) H , with T x = WXW^ and R A = [np(E)(E - A)]" 1 . Here U 
is the unitary matrix of eigenvectors of H and A = diag{Ai, . . . , Atv} stands for the 
diagonal matrix of the corresponding eigenvalues. The averaging over H then can be 
performed in two steps, the first step being the averaging over the Haar measure on 
the unitary group U(N). As this is again a special case of the IZHC integral it can 
be done explicitly. The important new feature however is that the N x N matrix T x 
is of a reduced rank, with its M <^C N nonzero eigenvalues coinciding with eigenvalues 
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7X C , c = 1, . . . , M of the matrix XW^W = jX, the rest of N — M eigenvalues being 
exactly zero. At the same time the resolvent matrix R\ is of the full rank N. The 
problem of performing the IZHC integral for two matrices of different rank can be most 
efficiently done by employing equation (A4) of the Appendix A in the paper [23] (which 
is in fact closely related to the so-called duality IZHC relation, see equation (17.3.8) in 
[22]). In our case it takes the form: 

<e X p H T, WW)), oc ^^ I A { Y) n j^^ *. ■ ■ ■ V* (*) 

where the integration goes over the complex variables y\, . . . , yu along contours parallel 
to the real axis such that sgn(Imy c ) = — sgn(x c ). Now we should perform the next 
step of the ensemble average over the eigenvalues A of H entering via the resolvent R\. 
After rescaling y c — > Ny c and a simple rearranging in the integrand we can see that the 
eigenvalue- averaged right-hand side of ( )28l) is proportional to 

dGtXM ~ N [ A{y}e-^^^)(n det(E-A) \ dyi _ m (29) 



A i X ) Jc T V=T det (E 



nNp(E)y, ^ 



In the limit A^ — > cxd the integrals over y c can be straightforwardly evaluated by the 
saddle-point method, with the saddle-point values being given by yc = -£-■ This 
is justified as equation ( FlOl) ensures that the expected value in the integrand tends for 
A^ — > oo to a well-defined limit of the order of unity along contours in the vicinity of 
the chosen saddle point. Moreover, one can show that the saddle-point can be reached 
by deforming the original contours without crossing any singularities of the integrand. 
Furthermore, A{Y^ S ^} oc A{X} detX~( M ~ 1 > and the Gaussian fluctuations around the 
saddle-point value yield the factor det X N ~~ X . Taking all these facts together we see that 
( 129]) indeed reproduces the expression for J r /3=2 (X) from (jSJ), and hence in the fixed 
amplitude model (HI) the i^-matrix in the limit N — > oo has the Cauchy distribution 
with density ([3]). This result has an interesting corollary. If w c are chosen to be the first 
M columns of the N x N identity matrix, then W\E — H)~ 1 W is nothing else as the 
M x M block of the resolvent (E — H)~ x . Therefore for invariant ensembles of Hermitian 
random matrices H, finite blocks of the resolvent of H are Cauchy-distributed in the 
limit of large matrix dimension. 
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